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Quantum information theory (QIT) emerged in physics as standard technique to extract relevant 
information from quantum systems. It has already contributed to the development of novel fields 
like quantum computing, quantum cryptography, and quantum complexity. This arises the question 
what information is stored according to QIT in molecules which are inherently quantum systems as 
well. Rigorous analysis of the central quantities of QIT on systematic series of molecules offered the 
introduction of the concept of chemical bond and aromaticity directly from physical principles and 
notions. We identify covalent bond, donor-acceptor dative bond, multiple bond, charge-shift bond, 
and aromaticity indicating unified picture of fundamental chemical models from ab initio. 


Extension of information theory * for quantum systems, 
called quantum information theory (QIT)^”"^, leads to 
emergence of non classical correlations, entanglement.^ 
In the past decade various concepts of QIT have matured 
to widely used tools in quantum many body physics.^ 
Several well known quantities have been redefined in term 
of entanglement shedding new light in our understanding 
of quantum systems. 

In chemistry such concepts have appeared only 
recently^“^^ although correlations among components of 
a finite system, like orbitals, can be regarded as exchang¬ 
ing information among such parties. For example, single 
orbital entropy provides information about how much an 
orbital is entangled with the rest of the system. ' In addi¬ 
tion, two-orbital mutual information^determines how 
orbitals communicate with each other, i.e., it gives the 
correlation between two orbitals as they are embedded 
in the whole system. The mutual information, however, 
contains correlations of both classical and quantum ori¬ 
gin. Such central quantities to describe orbital correla¬ 
tions have already been applied recently by some of us to 
analyze copper-oxide clusters^^, to dissect electron cor¬ 
relation effects in bond-forming and bond-breaking pro¬ 
cesses. In addition, entanglement structures have 

also been determined in photosystem lE^ and orbital en¬ 
tanglement analysis of the Ru-NO bond in a Ruthenium 
nitrosyl complex has also been carried out.^^ In all these 
works, correlations among the orbitals were measured 
in terms of quantum information entropies which are 
weighted averages of the eigenvalue spectrum of the var¬ 
ious subsystem density matrices (see Eq.(l)). Therefore, 
the more detailed information encoded in the the eigen¬ 
value spectrum^^’^^ and in the structure of the corre¬ 
sponding eigenstates of reduced density matrix has been 
lost. 

In this work, we present a systematic analysis of (multi- 
)orbital entanglement together with the probability dis¬ 
tribution of eigenstates given by the corresponding re¬ 
duced density matrices for all possible realizations of the 
two-orbital subsystems in series of handful of molecules. 
As a result, we find strong connection between our ap¬ 
proach and basic chemical models which allow us to de¬ 


scribe covalent bonds, donor-acceptor dative bonds, mul¬ 
tiple bonds, charge-shift bond, and aromaticity from a 
unified point of view. 

Theoretical background: When a system is split into 
two parts called a bipartite system (the two parts of¬ 
ten called Alice and Bob), the Hilbert space is H = 
If the system can be described by a pure 
state, the wave function is a linear combination of the 
tensor product of the basis functions of the two subsys- 

terns i-e, 1^-) = ® Up), 

where Cap is a complex matrix in general. The cor¬ 
relations between the two subsystems is fully quantum 
mechanical and called entanglement. The wave func¬ 
tion can also be written as a single sum of the prod¬ 
uct of transformed basis states due to Schmidt decom- 
position, i.e., = Ea=i ® where 

r < mm{dim'H^^\ dimH^^^), uja > 0, and = 1- 

The square of the \/^a Schmidt values are also equal 
to the eigenvalues of the so-called reduced density ma¬ 
trix, formed by tracing out one subsystem^, i.e., 

= TrB|v|/)(^-|, thus pfl = and 

are the eigenstates of and p^^\ respectively. 
If r = I the wave function |^) is a product state and 
the system is called separable, otherwise it is said to be 
entangled. 

In general, both subsystems are in a mixed state and 
the total system cannot be written as a single product 
of the states of the two subsystems. There are various 
quantities introduced to measure the mixedness of the 
subsystems and the strength of entanglement^ but they 
all must fulfill an important criterion namely entangle¬ 
ment monotonicity^^”^^. This means that the quantity 
cannot increase by local operations and classical commu¬ 
nication (LOCC). A widely used quantity is the entan¬ 
glement entropy given by the von-Neumann entropy of 
the reduced density matrix calculated as 

S{p^^^) = ( 1 ) 

a 

There are various possibilities to split a system into 
two or several subsystems. If subsystem (A) contains 


2 


a single orbital and subsystem (B) the rest of the or¬ 
bitals, the single-orbital entropy, can be cal¬ 
culated with i = where N is the number of 

orbitals/ Assuming four basis states per orbital, i.e., 
G {|0), |i). It), fhe theoretical maximum of 

is In4 with uja^ = 1/4 for all a = 1,...,4. 

The correlation between two orbitals, i and j, as they 
are embedded in the whole system is given by the two- 
orbital mutual information^ 

/bd) = S'(p(d) + 5'(p(i)) _ , ( 2 ) 

where are the single-orbital and 

two-orbital entropies, respectively. In this case is 

also a mixed state thus the mutual information contains 
correlations of both classical and quantum origin.An 
eigenvalue of denoted by corresponds to eigen- 
vector = T,ai,aj Cauajio^Mai) ® \<Paj) where ai 

and aj run from 1 to 4 and a = 1,... 16. Again the 
theoretical maximum of = In 16 corresponds to 

maximally entangled pure two-orbital state with = 
[1, 0..., 0]. This also means that orbital i and j are in 
maximally mixed state with cja^ = uja^ = 1/4 for all a 
and 5'(p^*^) = = In4, but orbital pair state (i^j) 

is in a pure state with 5'(p^^’-^^) = 0. 

In case of quantum chemical systems finite number 
of electrons are correlated on finite number of orbitals 
and the number of electrons, n, as well as the total 
spin projection, are good quantum numbers, 
also commute with the Hamilton operator and it has a 
block diagonal structure ^^’26,27 eigenstate 

of the two-orbital subsystem with quantum number pair 
(n, s^) (with n = 0,..., 4 and G { — 1, —1/2, 0,1/2,1} 
can be written as a linear combinations of basis states 
in the corresponding subspace. The set of coeffi¬ 

cients corresponding to basis states with quantum num¬ 
ber pair (n,s^) will be labeled by a vector Ca'^\n,s^) 
in order to use a compact notation for the rest of 
the paper. Therefore, c2’'^^(0,0) corresponds to |0,0), 
cF(l,-i) to {|04), ir 0 )}, cF(l, i) to {| 0 ,t), |t, 0 )}, 

0 /^/ 2 , 0 ) to {|o,n), It,;), i;,t), it;,o)}, c/^/ 2 ,- 1 ) to 

i; 4 ),c/^^( 2 ,i) to |t,t),c/^^(3,-i)to{|t;,;), i;,t;)}, 

i) to it;,t), it,t;)}, and c/'^^(4,o) to it;,t;)- 

In the rest of the paper, eigenvalues will be or- 

dered decreasingly, i.e, the largest value will correspond 
to (a = 1 and the smallest to a = 16. We will use the 
term relevant eigenvalue if it is one or two orders of mag¬ 
nitude larger than the remaining ones. The elements of 
the Ca'^\n, s^) vector will be given in terms of its largest 
element in order to show the ratio among the coefficients. 
Furthermore, in some cases the eqality among the vari¬ 
ous Ca'^\n, s^) vectors discussed in the next section will 
hold up to a spin reversal. 

Results and Discussion 

Fig. 1 shows the mutual information results of eight small 
molecules, BH 3 , CH 4 , NH 3 , H 2 O, HF, F 2 , N 2 , and CO. 


Numbers on each pictures only refer to different sites 
which orbitals are also shown around. Mutual informa¬ 
tion is indicated by the color of the lines between the 
sites. The lack of strong visible line between sites is the 
sign of negligible communication; these orbital correla¬ 
tions do not contain relevant information. The maximum 
of the mutual information is marked by the theoretical 
limit. In 16 ^ 2.77, as it was discussed earlier. 

The mutual information analysis of BH 3 reveals three 
strong entanglements given by = /(^T) = /(4,8) _ 

2.68 (Fig. la). All other 25 possible connections are zero 
within numerical accuracy. Therefore, site #2 and #6 
do not communicate with other sites which is in accor¬ 
dance with the chemical insight as site 7^2 is the in¬ 
ert Is core orbital while site 7^6 is the empty 2 p or¬ 
bital of the boron. The calculated mutual informa¬ 
tion of the correlated sites are identical and very close 
to the theoretical upper limit as their two-orbital re¬ 
duced density matrix is dominated by a single eigenvalue 

0.992. In the corresponding 
eigenvector, elements | 0 ,ti), It, i), It, t), Itt, 0 ), have 
almost the same relevance in the description of the in¬ 
teraction, i.e., c^^’^^( 2 , 0 ) = c^^’^^( 2 , 0 ) = c^^’^^( 2 , 0 ) = 
0.541 X [0.80,-1,1,0.88]. 

Similar picture is found for CH 4 molecule. Core or¬ 
bital, site 7 ^ 3 , does not entangle with other sites (Fig. lb) 
while there are four pairs of sites with considerable mu¬ 
tual information = /(^T) = /(4,8) _ j(5,9) _ 

2.68. Orbital pictures on Fig. lb show that in all 
four cases one hydrogen and carbon orbitals present 
strong mutual information. The corresponding relevant 
eigenvector elements are similar to the results of BH 3 
c/®^(2,0) = cf’^^(2,0) = c/®^(2,0) = c/®^(2,0) = 
0.543 X [0.79,-1,1,0.87]. Mutual information analysis 
of other small molecules like NH 3 , H 2 O, HF, and F 2 
(Fig. Ic-f) shows similar overall picture to BH 3 and CH 4 
indicating 3, 2, 1, and 1 very strong mutual informa¬ 
tion between sites, respectively (see also in SI I). These 
correlations always connect orbitals on different atoms 
and thus can be assigned to chemical bonds. The cal¬ 
culated mutual information of NH 3 , H 2 O, HF, and F 2 
is, however, only 2.52, 2.23, 1.85, and 1.47 smaller for 
these correlations than that in the case of BH 3 and CH 4 . 
Deviation from the theoretical limit reflects the relative 
difference of the elements of the eigenvector and the pres¬ 
ence of additional correlations of sites. For NH 3 , one 
eigenvalue dominates = 0.967 

with moderate difference in the weight of the eigenvec¬ 
tor elements, c^^’^^( 2 , 0 ) = cp’^^( 2 , 0 ) = c^^’^\ 2 , 0 ) = 
0.541 X [0.73,-1,1,0.94]. 

For H 2 O, the eigenvector elements have even 
more asymmetric distribution compared to NH 3 

(c/®^ (2,0) = (2,0) = 0.539 x [0.69, -1,1,0.98]). 

Besides, site 7 ^ 4 , the lone pair in the plane, has also 
minor contribution to the mutual information picture, 
/(4-6) = J(4,7) = 0.30 and = /(4,5) = o.22. These 

values are one order of magnitude smaller than the 
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FIG. 1. Graphical representation of mutual information for basic small molecules: BH 3 (a), GH 4 (b), NH 3 (c), H 2 O (d), HF 
(e), F 2 (f), N 2 (g), and GO (h). The color of the bonds according to the color scale (i) indicates the mutual information of the 
connected two sites. The shape of the molecular orbitals (isovalue 0.05) are shown directly near the sites. The blue and red 
colors of the orbitals refer to the sign of the wave function. Atoms are labeled with their usual colors; H, B, G, N, O, and F 
atoms marked with white, brown, tan, blue, red, and purple colors, respectively. 


usual ones = 2.23. Further analysis 

reveals that this secondary effect comes from the three 
electron regime given by the second largest eigenvalue, 
= 0.249 and the corresponding vector elements of 
^(4,6)(3 1^1)^ These secondary correlations are present 
between a site with two electrons and the two sites 
with large mutual information which may be regarded 
as hyperconjugative effect between a lone pair and an 
adjacent bond. Similar effects are observed for HF 
and F 2 as well which is also partly responsible for 
the lower value of I which is assigned to the chemical 
bond. However, the main reason is the large value of 
the second largest eigenvalue. For HF, 

= 0.118, while for F 2 = 0.174. It 


turns out that this second and third eigenvalues are 
associated with the three-electron regime. For HF, 

= 0.882 X [l,-0.53], for Fa 
^( 5 , 10 )^ 2 ^ ±i) = ± 5 ) = 0.707 X [1,1]. It formally 

means that an extra electron resonates between the two 
sites which is consistent with picture of charge-shift 
bond.^^^^ The orbital pictures of Fig. 1 underline this 
correspondence. 

Multiple bonds: We investigate N 2 and CO as model 
compounds for multiple bond systems (Fig. Ig, h). In¬ 
deed, the mutual information analysis results in three 
bonds for both molecules as expected but detailed in¬ 
vestigation identify fundamental differences. For tt- 
bond in N 2 , the eigenvector elements, c^^’^^( 2 , 0 ) = 
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FIG. 2 . Graphical representation of mutual information for the series of NH 3 BH 3 (a), NH 2 BH 2 (b), NHBH (c), GH 3 GH 3 (d), 
GH 2 GH 2 (e), GHGH (f), SiH 3 GH 3 (g), SiH 2 GH 2 (h), and SiHGH (i). The color of the bonds indicates the mutual information 
of the connected two sites. The shape of the important molecular orbitals (isovalue 0.05) are shown on the two sides where 
numbers near the orbitals indicate the site number. The blue and red colors of the orbitals refer to the sign of the wave function. 
Atoms are labeled with their usual colors; H, B, G, N, and Si atoms marked with white, brown, tan, blue, and yellow colors, 
respectively. 


^(3,10) ( 2 ^ 0) = 0.573 X [0.72,1, —1, 0.72], are similar to the 
previously discussed covalent bonds with the dominance 
of electron-sharing components. While the tt bonds of 
CO, however, have strong asymmetric distribution in 
eigenvector elements, c^^’^\2,0) = c^^’^\2,0) = 0.644 x 
[0.36, 0.80, —0.80,1], which surpass the usual asymmetry 
resulted from the polarization of the bond. Therefore, we 
conclude that QIT analysis may be able to differentiate 
between covalent and donor-acceptor bonds as well. 

Dative and multiple bonds: Unfortunately, owing to 
the symmetry of the CO molecule we cannot exam¬ 
ine separately the two tt bonds (one covalent and one 
donor-acceptor bonds) and support our hypothesis. To 
gain deeper insight, we investigate a series of molecules 
which can help to elucidate this question. We choose 
NH 3 BH 3 , NH 2 BH 2 , and NHBH for donor-acceptor test 
systems; CH 3 CH 3 , CH 2 CH 2 , and CHCH for isoelectronic 
apolar covalent reference systems. We also calculate 
SiH 3 CH 3 , SiH 2 CH 2 , and SiHCH as polarized covalent 
analogs which also serve as an example for bonds con¬ 
taining heavier element (see Fig. 2). In the rest of the pa¬ 
per, we only focus on the newly emerged bonding modes, 
while the previously discussed bond types are considered 


as known, if it is not mentioned otherwise. As an exam¬ 
ple, in this section, all C-H, N-H, and B-H bonds show 
similar picture as we have discussed for these type of 
bonds in the case of BH 3 , CH 4 , and NH 3 . We note that 
this result also indicates that QIT based analysis is ro¬ 
bust, yields the same relevant information for the same 
chemical moiety. We provide all details for all investi¬ 
gated molecules in Supporting Information H. 

QIT analysis for NH 3 BH 3 (Fig. 2a) provides the sim¬ 
ilar results to the donor-acceptor bond of CO. The mu¬ 
tual information of the donor-acceptor bond is signifi¬ 
cantly lower than for the covalent bond, /Ud3) _ 1 75 ^ 
The eigenvector elements show drastic asymmetric dis¬ 
tribution and the weight of electron-sharing parts is also 
smaller than in the previous examples: c^^’^^^( 2 , 0 ) = 
0.762 X [0.23,-0.58,0.58,1]. For the isoelectronic ana¬ 
log CH 3 CH 3 (Fig. 2 d), the C-C bond, site #5 and #13, 
has similar values to other covalent bonds; /U43) _ 
2.49, the eigenvector elements indicate difference from 
donor-acceptor bonding mode c^^’^^^(2,0) = 0.542 x 
[0.84,1,-1,0.84]. SiH 3 CH 3 with polarized Si-C cova¬ 
lent bond shows the expected results (Fig. 2g); mutual 
information is larger _ 2.28) than that of the 
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FIG. 3. Graphical representation of mutual information for the series of delocalized systems: GH 2 GHBH 2 (a), GH 2 GHNH 2 
(b), butadiene (c), cyclobutadiene (d), borole (e), benzene (f), pyrrole (g), furan (h), and thiophene (i). The color of the bonds 
indicates the mutual information of the connected two sites. The shape of the important molecular orbitals (isovalue 0.05) are 
shown on the two sides where numbers near the orbitals indicate the site number. The blue and red colors of the orbitals refer 
to the sign of the wave function. Atoms are labeled with their usual colors; H, B, G, N, O, and S atoms marked with white, 
brown, tan, blue, red, and yellow colors, respectively. 


donor-acceptor analog but somewhat smaller than that 
of the homonuclear analog. The important eigenvector 
elements show the expected asymmetry = 

0.553 X [1, —0.96, 0.96, 0.64]. Such direct comparison can 
be flawed by the large electronegativity difference of 
atoms therefore we investigate NH 2 BH 2 . The mutual 
information analysis suggests two strong interactions be¬ 
tween N and B atoms (Fig. 2b), and 

accordance with the double bond structure. Interest¬ 
ingly, the QIT results of these bonds are quite differ¬ 
ent. On one hand, 2.38 close to the values of 

previously mentioned covalent bonds and the eigenvec¬ 
tor elements of the largest also supports the po¬ 

larized covalent bond assignation based on previous ex- 
amples cf’'‘^(2,0) = 0.628 x [1,0.80,-0.80,0.51], On the 
other hand, /O.n) = 1,94 jg much lower and the eigen- 
vector elements is determined by the strongly asymmet- 
ric distribution cf’^^^(2,0) = 0.768 x [0.16,0.58, -0.58,1] 
similar to the donor-acceptor dative bond in NH 3 BH 3 . 


Assignation based on orbital images proves that the 
aforementioned results are consistent with the chemi¬ 
cal picture; the tt bond is the dative bond while the a 
is the covalent bond in NH 2 BH 2 . Although, these re¬ 
sults indicate no clear cut between covalent and dative 
bonds but the difference is prominent which is enough 
to assign them even within a double bond. Analog 
molecules CH 2 CH 2 (Fig. 2e) and SiH 2 CH 2 (Fig. 2h) 
show the expected results indicating covalent bonds; 
Ii^d) for these four bonds is in the range of 2 .3-2.5 
(see details in Supporting Information HE, II H). For 
NHBH (Fig. 2c), a sigma bond is given by = 

2.48, cf’^^(2,0) = 0.602 x [1,0.84,-0.84,0.60], while the 
two TT bonds are identical, because of symmetry rea- 
sons (J(5-9) = J( 6 . 10 ) = i.9r cf’®^( 2 , 0 ) = = 

0.644 X [0.35,0.80,-0.80,1]). These results are nearly 
the same as the results of CO suggesting the same bond¬ 
ing mode. For analog CHCH (Fig. 2 f), we obtain al¬ 
most the same results as for N 2 ; for the a bond: _ 
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2.61, cf’^^(2,0) = 0.516 X [0.94,1, -1,0.94], and for the tt 
bonds: J(®h 0 ) = 7 ( 6 , 9 ) =2.03, (2,0) = cf’®^(2,0) = 

0.571 X [0.73, —1,1, 0.73]. The same is valid for the anal¬ 
ysis of SiHCH (Fig. 2i) with the expected slight polariza¬ 
tion (see in SI III). 

Delocalized systems: After the successful determina¬ 
tion of multiple and dative bonding modes, we extend 
our investigation toward delocalized systems. As a tran¬ 
sition, we continue with the examination of CH 2 CH 2 BH 2 
and CH 2 CH 2 NH 2 molecules (Fig. 3a and 3b, respec¬ 
tively). Mutual information analysis of CH 2 CH 2 BH 2 re¬ 
veals a cyclic structure between site #2, #16, and #19. 
There is strong entanglement between site #2 and #16, 
7 ( 2 , 16 ) ^ 2.06, cf’^®^(2,0) = 0.594 x [0.69,-1,1,0.61], 
which corresponds to the C-C tt bond. While the other 
two mutual information data are an order of magnitude 
smaller, /(^d9) = 0.37, = 0.43. Further analy¬ 

sis points out that the main component arise from the 
one electron regime: (1, ±^) = 0.947 x [—1,0.34], 

= 0.971 X [1,0.25]. This is in accordance 
with the conjugated picture; the vacant orbital of the 
boron interacts with the pT^-orbital carbon atoms and 
forms a 3-center, 2-electron bond. Similar effects are 
found for CH 2 CH 2 NH 2 but the secondary interaction 
stem from the three electron regime: cp’^^^(3,±|) = 
0.934 X [—1,0.38] (see details in SI IIIB) as expected for 
a 3-center, 4-electron bond. 

To investigate longer delocalized systems, we choose 
butadiene (Fig. 3c). There are four interconnected sites 
in the mutual information picture. Two strong correla¬ 
tions, /<^^d9) _ j(i3,23) _ 2.13, are associated with the 
C-C TT bonds and there are secondary effects between 
them given by /(19,23) _ q /(3,i3) _ q which is 

consistent with the chemical picture of the delocalized tt- 
system. Analyzing longer delocalized molecule like hex- 
atriene shows similar results; there are secondary effects 
between the neighboring strong C-C tt bonds (see in SI 
HIE). 

Aromaticity: Molecules with cyclic delocalization have 
special place in chemists’ mind as the subject of the 
concept of aromaticity and antiaromaticity. Therefore, 
we extended our investigation to this direction. From 
butadiene and hexatriene one can derive the prototype 
antiaromatic and aromatic systems, cyclobutadiene and 
benzene (Fig. 3c and 3e, respectively, and in SI IV 
and V). For cyclobutadiene, mutual information anal¬ 
ysis indicates two strong correlations in the 7r-system, 
j(2,i8) _ j(i 2 , 22 ) _ 2 . 33 . Interestingly, the secondary 
effects observed in butadiene is disappeared in cyclobu¬ 
tadiene. We have not found any communication between 
the TT-subsystems within numerical accuracy, only the a- 
system show some minor effects probably due to their 
strained structure. Another antiaromatic compound bo- 
role shows the similar effects (Fig. 3e). In the 7r-system 
two C-C-B moieties are found, similar to CH 2 CH 2 BH 2 . 
However, there is no secondary effect in the butadiene 
moiety. Opposite effects can be seen for benzene. The 


strong mutual information between C-C 7 r-bonds become 
less dominant while secondary effects are even stronger 
compared to hexatriene and thus a cyclic structure is seen 

in the mutual information picture with very low value 
j(2,9) _ j(9,30) _ j(30,35) _ j(35,15) _ j(15,26) _ j(26,2) _ 

0.892 ± 0.007, (Fig. 3f). Interestingly, other weak inter¬ 
actions are found for the opposite sites, at para position 
of the aromatic ring, = /(9d5) _ j(26,30) _ q 

indicating their direct relationship. We have also investi¬ 
gated the aromatic series of furan, pyrrole, and thiophene 
(Fig. 3g-i). Interestingly, five-membered cyclic structure 
dominates the mutual information picture of the pyrrole 
TT-system which is in wide contrast to borole. The sec¬ 
ondary effects increased between butadiene moiety while 
the CH 2 CH 2 NH 2 structure is less emphasized. Cyclic de- 
localization with reduced mutual information between C- 
C double bonds and increased secondary correlations are 
found for furan and thiophene as well indicating similar 
effects to benzene and pyrrol. In the mutual informa¬ 
tion picture of furan and thiophene an additional weak 
structure is also found which are the hyperconjugative 
interaction of the lone pair with the adjacent cr-bonds 
similar to what has been found for water. 

In conclusion, we have introduced a novel approach to 
extract information from molecules based on QIT anal¬ 
ysis. Systematic investigation of handful of molecules 
using localized orbitals offers the introduction of the 
concept of chemical bond and aromaticity. We have 
shown on several examples how different chemical models 
like covalent bond, donor-acceptor dative bond, multiple 
bond, charge-shift bond, conjugation, and aromaticity 
follows from QIT. The discussed results indicate the 
unified picture of chemical concepts and therefore 
can help to elucidate their fundamental features and 
may lead to an improved definition of chemical bond. 
This study also closes the gap between state-of-the-art 
physical and traditional chemical models showing their 
mutual origin. Although there are many open questions 
regarding the QIT analysis of molecules we envision that 
our approach can be used in the future alternatively 
or together with well known chemical bond analysis 
methods^°“^^ and elucidate unique bonding modes. 

Methods 

In this paper we have performed numerical calcula¬ 
tions using the quantum chemistry version of the density 
matrix renormalization group (QC-DMRG) method. 

We have controlled the numerical accuracy using the 
dynamic block state selection (BBSS) approach^^ and 
the maximum number of block states varied in the range 
of 500-2000 for an a priory set quantum information loss 
threshold value y = 10“^. The ordering of molecular or¬ 
bitals along the one-dimensional topology of the DMRG 
was optimized using the Fiedler approach^^d^ 
active space was extended dynamically based on the 
dynamically extended active space (DBAS) procedure.^ 

Geometries have been optimized at HF/STO-3G level 
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of theory which yielded sufficient geometries in accor¬ 
dance with higher level methods. In QIT, we are es¬ 
pecially interested in the communication of subsystems 
which may communicate through active sites belonging 
to different subsystems. In molecules, atoms seem a nat¬ 
ural choice for the definition of subsystem. To create 
sites which correspond to one atom, we have applied lo¬ 
calized orbitals as sites. We have chosen Pipek-Mezey 
procedure^^ with tight threshold 10“^^ and minimized 
the number of atomic orbitals contributed in each local¬ 
ized orbitals. Therefore, they can be easily identified as 
part of a subsystem and because of the minimal basis set 
their chemical meaning is also clear which we can use for 
later comparison. All localized orbitals have been used 
in the DMRG procedure thus, as a result, we have car¬ 


ried out calculations at the FCI limit for all molecules. 
Then, results at the FCI limit have been analyzed in the 
paper. All preliminary calculations have been done by 
MOLPRO Version 2010.1^^. 
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In the Supporting Information (SI) we present the graphical representation of the single-site en¬ 
tropy and the correspondig cja^ eigenvalue spectrum (a = 1... 4) as a function of orbital 

index (i). Furthermore, we show pair-wise elements of the two-orbital mutual information 
indicated by the color of the lines between the sites. The lack of strong visible line between sites is 
the sign of negligible communication; these orbital correlations do not contain relevant information. 
The maximum of the mutual information is marked by the theoretical limit, In 16 2.77. The 

largest >0.1 eigenvalues of the corresponding two-orbital reduced matrix are collected in 
a table format together with the eigenvector coefficients obtained in the different quantum number 
sectors of the two-site subsystem. 
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II. DATIVE SYSTEMS 
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B. NH2BH2 
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C. NHBH 
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D. CH3CH3 
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E. CH2CH2 
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G. SiHaCHa 
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III. DELOCALIZED SYSTEMS 


A. CH2CH2BH2 



Orbital index 


i 

j 

I(iJ) 

Uoi >0.1 

II 

to 

Cfi 

II 

0 


t 

t , i 

ti ,- 

11 

20 

2.633 

0.985 

- 0.427 

- 0.541 

0.541 

- 0.483 

9 

10 

2.625 

0.984 

0.426 

- 0.541 

0.541 

0.483 

1 

12 

2.538 

0.968 

- 0.432 

- 0.543 

0.543 

- 0.474 

13 

14 

2.532 

0.967 

- 0.476 

- 0.542 

0.542 

- 0.430 

7 

15 

2.491 

0.958 

- 0.435 

- 0.543 

0.543 

- 0.471 

8 

17 

2.459 

0.955 

- 0.397 

0.529 

- 0.529 

- 0.532 

4 

5 

2.396 

0.933 

- 0.466 

0.524 

- 0.524 

- 0.484 

2 

16 

2.064 

0.866 

- 0.409 

0.594 

- 0.594 

- 0.359 

16 

19 

0.434 

0.189 

0.063 

- 0.321 

0.321 

0.889 

2 

19 

0.367 

0.187 

- 0.062 

- 0.437 

0.437 

- 0.784 


i 

j 

I(iJ) 

UJa >0.1 

n = 0 

n = 3 , s * = - 1/2 

n = 3 , = 1/2 

n = 4 : 


tl , i 


ti , t 

t , 

ti , 

16 

19 

0.434 

0.116 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

2 

19 

0.367 

0.150 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 


OI8 


20 


11 




6 

O 



19 ° <\ 1 ^ 0-^15 

\ ^ 

2 12 


i 

j 

I(iJ) 

UJa >0.1 

n = l , = 1/2 

n = 1, = -1/2 

t ,- 

t 

i ,- 


16 

19 

0.434 

0.343 

0.343 

- 0.947 

0.000 

0.321 

0.000 

0.000 

0.947 

0.000 

- 0.321 

2 

19 

0.367 

0.327 

0.327 

0.971 

0.000 

0.240 

0.000 

0.000 

- 0.971 

0.000 

- 0.240 





















































































19 
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D. Cyclopentadiene 
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E. Hexatriene 
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IV. ANTIROMATICITY 


A. Cyclobutadiene 
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B. Borole 
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V. AROMATICITY 
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B. Pyrrole 
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